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Abstract 

We present a first numerical study of lattice QCD with 0(a) improved Wilson 
quarks and a chirally twisted mass term. Renormalized correlation functions are de- 
rived from the Schrodinger functional and evaluated in an intermediate space-time vol- 
ume of size 0.75 3 x 1.5 fm 4 . In the quenched approximation precise results are then 
obtained with a moderate computational effort, allowing for a detailed study of the 
continuum approach. The latter is discussed in terms of observables which converge to 
meson masses and decay constants in the limit of large space-time volume. In the 0(a) 
improved theory we find residual cutoff effects to be at the level of a few percents at 
a ~ 0.1 fm. 
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1 Introduction 



Lattice twisted mass QCD (lattice tmQCD) has been introduced in Refs. [T],^] as a so- 
lution to the problem of spurious quark zero modes, which plague lattice computations 
with light quarks of the Wilson type, especially if the action is 0(a) improved. The 
occurrence of spurious quark zero modes causes a breakdown of the quenched and par- 
tially quenched approximations, as well as technical problems in the fully unquenched 
simulations. In Ref. || it has been shown how Symanzik's on-shell improvement pro- 
gramme can be implemented in the framework of this new lattice regularization, 
which is intended for QCD with two mass degenerate quarks. 

The purpose of this paper is to investigate the scaling violations in lattice twisted 
mass QCD, following the lines of the study presented in Ref. || for the 0(a) improved 
Wilson lattice action. The tmQCD lattice regularization differs from the latter only 
by the parameterization of the quark mass term, which is rotated in the chiral flavour 
space, whereas the Wilson term remains in the standard form. This difference is the 
key to avoid spurious quark zero modes and can viewed as a change of the quark field 
basis that leaves unchanged -up to cutoff effects- the physical content of the theory. 
Therefore it often happens that a certain physical quantity is obtained from a different 
correlation function -with different cutoff effects- compared to the case of standard 
quark mass parameterization. Moreover all the cutoff effects that are proportional to 
(some power of) the quark mass may quantitatively change, although for light quarks 
this is expected to be a small effect. These remarks motivate our investigation of scaling 
violations in lattice tmQCD. 

The physical parameters of the present study have been chosen so as to be in a 
situation similar to that of the scaling test of Ref. Namely, we consider a system 
of finite size, (L 3 x T ~ 0.75 3 x 1.5) fm 4 , with Schrodinger functional boundary con- 
ditions, and give the relevant renormalized quark mass parameter, Mr, a value such 
that LMr ~ 0.15. For such a system, we study the approach to the continuum limit 
of a few renormalized observables, which in the limit T — > oo have the same physical 
interpretation as the observables studied in Ref. ||. When also the limit of large L 
is taken, the observables turn into the pion mass, the p-meson mass, the pion decay 
constant and a quantity related to the p- meson decay constant. Based on the study of 
Ref. H in large spatial volume, we expect the cutoff effects observed at L = 0.75 fm to 
be indicative of the size of the lattice artifacts in infinite volume. 

In Section 2 we introduce the relevant correlation functions within the Schrodinger 
functional setup for lattice tmQCD. The renormalized observables of interest are con- 
structed in Section 3, where the renormalization scheme adopted for tmQCD is also 
specified. Section 4 presents numerical details and a discussion of our results, while 
conclusions are drawn in Section 5. A preliminary report on the present work has al- 
ready appeared in Ref. 0. In the following we assume that the reader is familiar with 
Refs. [|l],||]3| and refer to the equations of Ref. Q by using the prefix I. 
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2 Schrodinger functional correlation functions 



The Schrodinger functional (SF) for lattice tmQCD has been introduced in Ref. ||, 
where it is defined as the integral kernel of the integer power T/a of the transfer matrix. 
It admits the following representation: 

Z[p',p',C';p,p,C] = [ D[C7]DMD[fl e~ s ^\ (2.1) 



where S[U, vp, tp] is the Euclidean action of tmQCD and the arguments of Z are the 
prescribed boundary values of the gauge and quark fields at xo = (C, p, p) and 
xq = T (C, p', p'). As usual, the Dirichlet time-boundary conditions for the quark 
fields involve the projectors P± = (1 ± 7o)/2. 

Renormalizability and 0(a) improvement of the SF for tmQCD have been discussed 
in Section 3 of Ref. |J. We take over the outcome of that discussion and extend the 
action S[U,^, ip] to include all the counterterms that are needed for renormalization and 
0(a) improvement, as detailed in eq. (1.3.4). In particular, adopting the same notational 
conventions as in Refs. HJ8|, the quark action Sf[U, "tp, tp] takes the same form as on the 
infinite lattice: 

S F [U, iP] = a 4 Y, {D + 5D + m + ^ q75 r 3 ) ^(x), (2.2) 

X 

where 5D stands for the sum of the volume and the boundary 0(a) counterterms. 
The SF correlation functions can be written in the form 

(?) = Sz- 1 f D[E7p[$D[$ T e - s[u '^ ] ] , (2.3) 

I J J p'=p'=p=p=0;C'=C=0 

where T stands for any product of fields localized both in the interior of the SF box and 
on its time-boundaries. For instance, quark and antiquark fields at x$ = are given by 

C(x) = P_C(x) 



C(x)=C(x)P+ = -j^y (2-4) 
The reader is referred to [^J8| for any undefined conventions. 

2.1 Bare SF correlation functions 

Within this SF setup we now define a few on-shell correlation functions that involve 
quark-antiquark pairs of boundary fields and the following isovector composite fields: 



V«(x) = ^x) 7 4r>(x), 
P a {x) = ^(x) 75 irXx), 



T^(x) = ^(x)ia^T a ^(x) , (2.5) 

where a^ v = (i/2)[7 M , 7„]. 
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In addition to the /-correlators that were introduced in Ref. ||], 

ff(x ) = -{A a (x)O b ), 
ff(x ) = -{P a (x)O b ), 

f$>(x ) = -(V a (x)O b ), (2.6) 

we also consider some further correlation functions: 

3 



k=l 



k=l 
3 



k=l 



The isospin indices a and 6 are restricted to take values in the set {1, 2} for reasons to 
be explained below, while the boundary fields O a and are defined by: 



O a = a 6 ^TC(y)75^C(z), 

Qt = a^~C(y)lk\r\^) . (2. 



For the purpose of boundary field renormalization, we also need to evaluate a 
boundary-to-boundary correlator: 

ff = -^(0' a O b ), (2.9) 

where ,a is defined analogously to O a but with derivatives with respect to quark 
boundary fields at xq = T rather than at xq = 0. 



2.2 Flavour structure of the SF correlators 

As long as the isospin indices a and b take the values 1 or 2, one can show H that the 
lattice symmetries of the tmQCD SF imply some exact properties of the bare correlation 
functions. For the /-correlators these relations are summarized by eqs. (1. 3. 48)^(1.3.49), 
while for the correlator ff 1 one finds 

/i n = /i 22 , /i 12 = /i 21 = 0. (2.10) 

The analogous relations for the /c-correlators, which can be easily derived along the lines 
of Ref. §, read 

4 2 (^o) = 4 2 (x ) = k^(xo) = 0, (2.11) 

and for X = V, T, A 

Af (x ) = kg(x ), k$(x ) = -fc x 2 (x ) • (2.12) 
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The non- vanishing correlators to be evaluated in practice can hence be chosen as fjj(xo), 
fr(xo), /v- 2 (x ), k$(xo), kf{x ) and f\\ 

An explicit representation of the /-correlators in terms of the boundary-to-bulk 
quark propagators is given in eqs. (I.3.50)-(I.3.52). The analogous representations for 
the /c-correlators and f\ l readQ: 

ky(x ) = |(|^tr {H + {x)^ 5lk H + (x)w 5 }\ , 
\ k=l I G 

3 



kr(x ) = \ / \ ^2 tr |^+(x) t 75707fe^+(a;)7fc75} ) , 

\ k=l I g 

kf{x a ) = \l\Y. iT {H+(x)h k H + (xh k75 }\ (2.13) 



k=l 



and 



^2 „6 



y.z \ 



(2.14) 



y =zo=T-a i 



G 



where H+(x), eq. (1.3.43), is the first flavour component of the boundary-to-bulk quark 
propagator H{x) defined in eq. (1.3.38). We remark that ft 1 > 0. 

From the setup of the SF for tmQCD Q one can readily see that evaluating H + 
amounts to solving for < xq < T the one-flavour system 

(D + 5D + m + i W )H+(x) = cta^S^Uix - a6, 0)- l P+, (2.15) 

with the boundary conditions 

P + H + (x)\ xo=Q = P^H + (x)\ XQ=T = 0. (2.16) 

The solution H+(x) of eq. ( p.15 ) satisfies H + {x)P + = H + (x) and is trivially related to 
the boundary-to-bulk quark propagator H + (x): 

H + {x) = H + {x) - 8 X0fi P + . (2.17) 

Provided the triplet isospin indices are restricted to the values 1 and 2, we are able 
to express the correlation functions of Subsection 2.1 in terms of H + (x) alone, which 
saves about a factor of two in CPU-time. Since the full physical isospin symmetry 
is expected to be restored in the continuum limit of lattice tmQCD Q, the above 
restriction implies no loss of physical information. 



3 Renormalization scheme and scaling observables 

The SF for lattice tmQCD is expected to be ultraviolet finite after renormalization of 
the bare parameters in the action, <7q, mo and /U q , and the boundary quark fields Q. As 

x As in Ref. | p2[ , the bracket (. . .)g means an average over the gauge fields with the effective gauge 
action. In the quenched approximation the average is performed with the pure gauge action. 
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the latter renormalize multiplicatively and are set to zero, here we do not have to worry 
about their renormalization. In the following we specify our renormalization scheme 
for the parameters in the action and the correlation functions. Provided that all the 
relevant improvement coefficients are given their proper values, the mutual relations 
among renormalized parameters and observables are free from 0(a) cutoff effects. 



3.1 Renormalized parameters 

Since we work in the quenched approximation to QCD, it is convenient to renormalize 
the gauge coupling by eliminating g$ in favour of the hadronic length scale ro ~ 0.5 fm 
H and then express all the physical quantities in units of ro- In large physical volume 
the relation between f3 = 6/g$ and a/ro has been evaluated [|l(]] with a relative accuracy 
of about 0.5%. Since for the present scaling study we are working in an intermediate 
volume, where finite-size effects are non- negligible, we keep constant the ratio L/tq, 

— = T— 1 (/?) - = 1.49, (3.1) 
r |_ r oJ a 

while approaching the continuum limit. We choose T/a = 2L/a and values of L/a 
such that the values of (3 lie in the range 6 < (3 < 6.5, which is the one of interest for 
quenched lattice QCD with the Wilson plaquette action. 

As for the renormalization of the quark mass parameters, we require 

Lm R = 0.020 , 

L/z R = 0.153, (3.2) 

where m R and are the renormalized quark mass parameters introduced in eqs. (1.2.8)- 
(1.2.9). Following Ref. ||, the renormalized twisted mass parameter fi R is related to 
the bare masses /i q and m q = mo — m c via 

fj, R = 2^(1 + b^am q )fi q , (3.3) 

where is an improvement coefficient introduced in ||. The exact lattice PCVC 
relation (1.2.14) implies that we can set: 

Z„ = Zp 1 . (3.4) 

The renormalization constant of the isotriplet pseudoscalar density Zp is evaluated in 
the SF scheme at the momentum scale q = (1.436ro) -1 , using the results of Ref. [15]. 



The value of m H in the same scheme and at the same scale is computed from the 
renormalized PCAC relation as discussed later on. 

We recall from Ref. M that in renormalized tmQCD the "polar" quark mass 



M R = y/ml + m| (3.5) 
plays the role of the renormalized quark mass. The angle a, defined by 

tana = — ^ , (3-6) 

m R 

can be chosen arbitrarily and just determines the physical interpretation of the tmQCD 
correlation functions. The numerical values on the r.h.s. of eq. (3.2) yield LM n ~ 0.154, 



which is in the range of the scaling study ||, and a value of a that is far from zero, 
namely n/2 — a ~ 0.130. 
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3.2 Renormalized and O(a) improved correlators 



/a + CAadofp - a/xq b A / v (x ) 



The definition of renormalized and 0(a) improved SF correlators is guided by the form 
of the renormalized and 0(a) improved bulk fields, eqs. (1.2.10)— (1.2.12) and 

(T R )% = Z T (1 + b T am q )[T^ + c T a(d^ a - d u V£)] , (3.7) 

where denotes the symmetric lattice derivative in the direction of the unit vector ft. 
We hence define: 

[/"Mr = [ZJ1 +Knm tl tfZ x U +h x 
[/v 2 (^o)]h = 

[/^ )]r = 
[A^Mr = 

[A#(*o)]r = 



[Z c (l + bc_am q )] 2 Zy{l + byam q ) fy + a/i q by f A (^o) , 
[Z ? (l + 6 ? am q )] 2 Z P (l + b P am q ) fp-(x ) , 

[Z{(1 + b^am q )] 2 Zy(l + byam q ) ky + cvadofcr 1 — a/i q oy fc\ 2 (zo) , 



[Z f (1 + 6^am q )] Z T (l + frTam q ) 



CTddok- 



(xo) , 



[*a(*o)]r 
and 



[Z c (l + b^arriq)} Z A (1 + b A am q ) k A + a/i q b A ky (x ) 



(3. 



[ft 



Hi 



[Z c (l + 6 c am q )] 4 /i n 



(3.9) 

The improvement coefficients b A and by have been introduced in || ; whereas all the 
remaining improvement coefficients are the same as in lattice QCD with standard quark 
mass parameterization. We remark that the expressions for [/y 2 ] R and [fc A 2 ] R are inde- 
pendent of cy and ca, respectively, because of the translational invariance of the theory 
in the spatial directions. 



3.3 Pion and p-meson channel correlators 

The definition of the observables for this scaling test is inspired by the criterion of 
considering observables that in the limit of large T and large L turn into the pion and 
/9-meson mass and decay constant, except for the normalization of the p-meson decay 
constant which is not the physical one. The same criterion was followed in the scaling 
study of Ref. §. 

The first step in the construction of the meson observables is to build linear com- 
binations of the correlators in eq. ( |3.8| ) so to yield at time xq insertions of operators 
with the appropriate quantum numbers to create/annihilate a pion or a p- meson (or 
higher states in the same channels). According to the relation between renormalized 
correlation functions of QCD and tmQCD in infinite volume |2|, such operators can be 
written as follows: 

(A^)g(x) = co S a(A R y (x)+e 3ac sma(V R ) c (x), 
(P R ) a (x) = (P R ) a (x), 

(V£)%(x) = cos a(V R y k (x)+e 3ac sin a(A R )%(x) , 
CMo(x) = (« (x). (3.10) 
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We remark that the expression of local operators with given physical quantum numbers 
in the tmQCD quark basis does not depend on the choice of boundary conditions, and 
the results of Ref. Q are hence valid in the present context. The situation is different 
for the correlation function themselves, so that the SF correlators defined below cannot 
be directly compared -at least for finite extent T of the SF- with those computed at 
the same value of M R and a = 0. A detailed discussion of this point is deferred to a 
forthcoming publication Jll[ |, 

The correlators containing the operator insertions in eq. dUOl) with isospin index 
a = 1 are correspondingly given by 

[/1'( 2; o)]r = cosa[/A 1 (x )] R - sina[/v 2 (x )] R , 

[ky,(xo)}n = cos a[ky (x )] K - sinafA^ 2 (z )] R > 
[A#(x )] R = [A#(*o)]r. (3.11) 
With the SF-boundary fields O 1 or Q\ introduced in Subsection 2.1, one expects that 



in the limit of large xq and T — xq the correlators in eq. ( 3.11 ) are dominated by the 



"pion" and "p- meson" states in spatial volume L 3 , respectively. 



3.4 The observables of this scaling test 



We are now ready to define the scaling observables which we will focus on in the re- 
maining part of this paper. 

• Meson observables 

In terms of the above correlators, eq. ( |3.11| ), the estimators of the finite volume 
pion (PS) and p-meson (V) masses read: 



m PS 



my 



[/£]h 



m PS 



x =T/2 



V'J R 



A'JR 
11 



x =T/2 



R 



mv = ~iav~ ■ (3 ' 12) 

x =T/2 

It should be noted that at finite T the quantities mps and my need not coincide 
with mps and my, respectively, because they may receive contributions from states 
heavier than the finite volume pion or p-meson. 

The estimators of the finite volume pion and p-meson decay constants read 



xo=T/2 
x =T/2 ' 
xo=T/2 ' 
x =T/2 ' 



??PS = 


[/i 11 


-1/2 
R 




VPS = 


[/i 11 


-1/2 
R 




W = 


[/i 11 


-1/2 
R 


Cy [k$(x )] B 


rjy = 


[/i 11 


-1/2 
R 


Cy [k Y \(x )] R 



(3.13) 
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where the normalization constants Cpg and Cy are given by 



C PS = . , Cy = , . (3.14) 



L 3 m ps yL 3 m^ 

The constants Cps and Cy are defined analogously in terms of mps and fhy. The 
normalization constant Cps is chosen such that rjps — ► F n as T = 2L — > oo. In 
the same limit the quantity rjy, due to its unphysical normalization, does not 
approach the (inverse) decay constant of the /9-meson, but one may expect the 
cutoff effects to be similar to those of the properly normalized estimator of 1/F p . 
Analogous remarks hold for the quantities rjp$ and rjy, which differ from ryps and 
r/y only by their normalization. 

PCVC and PCAC quark masses 

The PCVC and PCAC operator relations of renormalized tmQCD, which follow 
from the flavour chiral Ward identities ||, imply corresponding relations among 
the renormalized SF correlators introduced above: 

da[ti 2 (x )]n = -2fi R [f^(x )] R (3.15) 

and 

<9o[/IVo)]r = 2m R [f P A (x )] R • (3.16) 

As a consequence of the improvement of the bulk action and relevant operators, at 
finite lattice spacing a the cutoff effects on these relations are 0(a 2 ), even without 
improving the SF-boundary action and fields. 

A way of checking the size of the residual cutoff effects in the PCVC relation, 
eq. (|3.15| ), is to consider the quantity 

''PCVC = , (3.17) 
A* 

where p, is the estimate of the bare current twisted mass obtained from the SF 
correlators: 

So/y 2 ^) 



= —Zy(l + byarriq) 



2/£ 1 (zo) 



(3.18) 

x =T/2 



In the definition of \x we have left out the improvement coefficients that are only 
perturbatively known^. Close to the continuum limit we expect: 

rpcvc = 1 - am^Z- l [bp + b fl + ZZyby] + 0(a 2 ) , (3.19) 



where Z = Z m Zp/ZA (see e.g. Ref. [12] for a recent non-perturbative estimate of 
Z as a function of g 2 ). In view of the very small values of am R that correspond to 
L/a > 8 and Lm R = 0.020, the cutoff effects on rpcvc should be dominated by 
the terms 0(a 2 ), and the sensitivity to the combination [bp + b^ + ZZyby] ~ 0(1) 
should hence be very small. Our data (see Section 4) confirm this expectation. 



2 However, a non-perturbative estimate of bp at (3 = 6 and (3 = 6.2 has been given in Ref. [ ]l3| . 
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The PC AC relation, eq. ( 3.16| ), is instead exploited to evaluate the renormalized 
standard mass, which is defined by 



d [fi l (xo)]n 



2[/ I !, 1 (^o)] R 



x =T/2 



(3.20) 



One could of course evaluate m R by eq. (1.2.5) and eq. (1.2.8), but then the absolute 
accuracy on Lm R is essentially limited by the accuracy on m c = to c (<7q). Turning the 
argument around, one can exploit the determination of ra R obtained from eq. ( p. 20 ) 
to estimate the critical quark mass m c , which is independent of a up to cutoff effects 
H|]. Within the tmQCD regularization, the estimate of m c can be performed for any (3 
with no singularities in the computation of quark propagators by working at very small 
values of m R and reasonable finite values of Moreover from Ref. j|] it can be argued 
that an 0(a) improved evaluation of m c in lattice tmQCD requires the non-perturbative 
knowledge of a certain combination of the improvement coefficients 6a and b m : 

b m - (ZZy)~ l b A . 



Work in this direction is currently in progress by the Tor Vergata APE group [14|. 



4 Numerical details and results 



The basic idea of any scaling test is to approach the continuum limit along a line 
in bare parameter space where all renormalized parameters are kept constant, which is 
achieved here by the renormalization conditions specified in Subsection 3.1. Under these 
conditions the renormalized and (almost) 0(a) improved observables that we introduced 
in Subsection 3.4 are expected to depend on a/L only and converge to a well-defined 
continuum limit as a/L — > with (almost) no scaling violations linear in a/L. 



4.1 Renormalization constants and improvement coefficients 



The renormalization conditions for the gauge coupling, eq. ( [3.1D , and for the two quark 
mass parameters, eq. fl3.2|) , where the scheme dependence arises only from Zp, are suf- 
ficient to render ultraviolet finite the observables introduced in Subsection 3.4. Indeed, 
no observables depend on the boundary field renormalization factor Z^(l + b^am q ), as 
well as on the product Z^(l + bTam q ). Concerning Z& and Zy, which are needed to re- 
move lattice artifacts that vanish more slowly than a as a — > 0, we employ the available 



non-perturbative estimates from Ref. [17]. Since all these renormalization constants are 
defined either at the chiral point of quenched QCD or in the Yang-Mills SU(3) theory, 
we have actually set up a non-perturbative quark mass independent scheme. 

In order to further reduce the scaling violations of our observables, we have to give 
proper values to all the relevant improvement coefficients. In the limit of large time 
extent T of the SF system, the O(a) improvement of the bulk action and operators is 
sufficient to improve the scaling observables defined in Subsection 3.4. As we actually 



work at T ~ 1.5 fm, this statement remains valid only for the quantity rpcvC; e Q- (3-17) 



The remaining scaling observables, which are defined at xq = T/2, are indeed not 
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completely dominated by the pion or the p-meson state and hence do still depend on 
the details of the SF-boundary fields. As a consequence, the 0(a) improvement of the 
SF boundary action and fields can not be neglected. 

Let us start with the coefficients that are relevant for the improvement of the 
massless theory. We employ the non-perturbative estimates of c sw , ca and cy computed 
in Refs. |l^Jl^], while setting ct to its one- loop value[] pOfl . As for the improvement of 
the SF-boundaries, the only coefficients relevant for our study are ct and c*t, which are 



both set to their one-loop values [21 ,22] . We discuss below the impact of reasonable 
changes of these values on our scaling observables. 

In view of the very small values of am R corresponding to Lm R = 0.020, the im- 
provement coefficients multiplying counterterms that are linear in am R need not be very 



precisely tuned. Nevertheless, we adopted non-perturbative estimates of 6a — 6p |L2|] 
and by J2Q] , and one-loop estimates for 6a an d 6p |J20|| , whereas 6t is not necessary at 
all for the improvement of our observables. It should be noted that, due to the quark 



mass renormalization conditions (^2), which entail a ~ 7r/2, the correlator [/ a }(xo)]r 
depends on 6a only via a term proportional to am R cos a <C 1. For the same reason, 
the ratio ^[fcy^R/f&y^R an d thus my are almost independent of 6a- 

Among the improvement coefficients that multiply counterterms of order a/d B [||, 
in this scaling test we only need to properly tune 6a, eq. (1.2.10), by, eq. (1.2.11), and 
6i, eq. (1.3.6). Moreover, since the set of improvement coefficients { 6i, 6 m , 6^, 6a, by} is 
actually redundant, as explained in Ref. ||, one of them can be arbitrarily prescribed. 
For practical reasons that are specific to this scaling test, we find it convenient to 
set 6i(oq) = 1 exactly. We can then completely forget about the corresponding O(a) 
boundary counterterms. Concerning 6a, by and 6 M , they are set to the one-loop values 
|| that follow from the above prescription of b\, i.e. 

6 A = 0.0213^, &v = 0.0053 5 g, h = -0.0440 5o 2 . (4.1) 

We remark that the counterterm proportional to 6^ is of order am R , while 6 m is not 
needed at all for this study. It remains to be checked a posteriori whether the residual 
0(afi R ) effects are significant in comparison to the higher order scaling violations and 
the statistical uncertainties of our observables. 



4.2 Simulation parameters and analysis of the raw data 

In order to check the size of the scaling violations in our observables and the rate of 
the approach to the continuum limit, we perform simulations with four different lattice 
resolutions in the range 2GeV < 1/a < 4GeV, while enforcing the renormalization 
conditions detailed above. Throughout the whole analysis we adopt the definitions of 
the renormalization constants and the improvement coefficients specified in Subsection 
4.1 and evaluate them at the values of (3 = 6/oq chosen for our simulations. 

An overview of the bare parameters, the corresponding renormalized parameters 
and the accumulated statistics is given in Table [l]. The values of L/tq and L/x R follow 
from our choice of bare parameters, and the quoted uncertainties stem from the statis- 



tical errors on a/ro [|H| and on Zp = Zp(q) [15], respectively. The uncertainty on Lm 



' R 



3 We recall that ct is only relevant for the improvement of mv and 77V- 
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reflects instead the statistical error on the SF correlators in eq. (ET2CD , which in turn 
is a combination of the statistical errors on the bare SF correlators, as evaluated via 
simulations of lattice tmQCD, and the known ratio Zp^jZ^. Our Monte Carlo simula- 



set 


L/a 




a/iq 


K 


L/r 




Lm R 


-^rneas 


A 


8 


6.0 


0.01 


0.134952 


1.490(6) 


0.1529(8) 


0.0228(23) 


7680 


B 


10 


6.14 


0.00794 


0.135614 


1.486(7) 


0.1530(8) 


0.0203(30) 


2880 


C 


12 


6.26 


0.00659 


0.135742 


1.495(7) 


0.1530(8) 


0.0201(23) 


3072 


D 


16 


6.47 


0.00493 


0.135611 


1.488(7) 


0.1529(8) 


0.0180(24) 


1680 



Table 1: The bare and renormalized parameters for our data sets and the number 
(N meaB ) of computed SF quark propagators on decorrelated gauge backgrounds. 



tions of quenched lattice tmQCD are performed on the APE100 parallel computers with 
32-256 nodes at INFN Milan and DESY Zeuthen. The parallelization of the program 
and the machine topology allow us to simulate several independent replica of the smaller 
systems (i.e. A,B and C in Table [l]) at the same time. The computational effort needed 
for the present scaling test amounts to about 75 Gflops x days. 

The gauge configurations are generated using a standard hybrid overrelaxation 
algorithm. The single iteration is defined by one heatbath step followed by iVoR = 
L/2a + 1 microcanonical reflection steps. The correlation functions are evaluated by 
averaging over sequential gauge field configurations separated by 50 iterations. For 
the computation of the quark propagators H + (x) entering our observables, we use the 
BiCGStab inversion algorithm with SSOR preconditioning [16]. As a stopping criterion 
for the inversion algorithm we require the square norm of the dynamical residue, as 
defined in Ref. |16], to be 10 13 times smaller than the square norm of the solution. For 
the finest lattice spacing considered, convergence is always reached within a number of 
BiCGStab iterations between 80 and 120. 

A binning analysis of our data shows that for all simulation points consecutive mea- 
surements of the correlation functions (in the above specified sense) can be effectively 
taken as statistically independent. As our observables are non-linear combination of 
the basic correlation functions, we adopt a single-elimination jackknife procedure for 
the evaluation of their statistical errors. 



4.3 Continuum limit extrapolations 

By the above analysis procedure we obtain the results for rpcvc shown in Figure [l] and 
the values of the meson observables quoted in Table [2|. 

The errors on our results for rpcvc already take into account the small statistical 
uncertainty^] on Zy, whereas the product byam q <C 1 is taken with no error. Inspec- 
tion of Figure Q immediately reveals that the residual cutoff effects of order am R -see 
eq. ( 3.1G| )— are completely negligible with respect to the very small statistical errors and 



higher order scaling violations. The very tiny mismatch between the values of Lm R and 



We neglect here the tiny systematic uncertainties associated to the determination of Zv- see Ref. J17 
for details. 
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Figure 1: Our results for rpcvc an d their continuum limit, obtained by fitting the four 
data points to a first order polynomial in (a/L) 2 . 



the nominal value 0.020 is hence completely negligible in this case. Moreover, due to its 
very definition, rpcvc is independent of Zp and cy. We hence conclude that for a ~ ir/2 
and the values of a// q that are relevant for (quenched) QCD in the chiral regime, the 
PCVC relation shows surprisingly small cutoff effects and is effectively O(a) improved 
once the proper value of c sw is employed. 

In the case of our meson observables, which may be quite sensitive to mismatches 
and uncertainties in the renormalization conditions as well as to uncancelled O(a) cut- 
off effects stemming from the SF-boundaries, we have performed a slightly more refined 
analysis before producing scaling plots and attempting continuum extrapolations. The 
starting point of this further analysis is represented by Table ^, which is directly ob- 
tained from the simulation data by considering all the renormalization constants and 
improvement coefficients with no error. The quoted errors arise just from the statistical 
fluctuations of the observables over the samples of gauge configurations produced at 
the bare parameters of Table |l[ The remaining uncertainties are taken into account as 
follows. 

• Uncertainties on the renormalization and improvement coefficients 

The observables mps and my have a very tiny dependence on Z^jZy which dis- 
appears in the limit a — > ir/2. Since a ~ 7r/2, the would-be decay constants of the 
pion and the p- meson are almost proportional to Zy and Za, respectively. The 
statistical uncertainties on Zy and Za 0], which are of about 0.01% and 1% of 
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set 


'"'Fb J -' 


Uly L 


'/Fn- L ' 


'IV 


A 


1 866(1 4"l 


Zj . V7 s7 O ^ 1(J } 


541 9(34) 


1 528(21 1 


B 


1 805(21) 

1 • UU'J 1 " J- ) 


2 652(271 


5570(561 


1565(341 


C 


1 831 (21 


2 646(28) 


5514/531 


1615(36") 


D 


1.825(27) 


2.648(35) 


0.5510(70) 


0.1601(46) 


set 


m PS L 


myL 




rjy 


A 


1.713(8) 


2.398(12) 


0.5654(25) 


0.1817(21) 


B 


1.667(11) 


2.337(19) 


0.5793(38) 


0.1889(34) 


C 


1.680(10) 


2.345(19) 


0.5751(36) 


0.1933(34) 


D 


1.659(13) 


2.323(25) 


0.5779(49) 


0.1949(47) 



Table 2: "Raw" results for our meson observables: the quoted errors arise just from 
statistical fluctuations over our samples of gauge configurations. 



the mean values, respectively, are added quadratically to the errors in Table 



The systematic uncertainties on Zy and Za [17] are shown separately as errors on 
the continuum limit extrapolations, see Table |j. 

The uncertainties on the improvement coefficients that are needed to subtract ef- 
fects of order a/x q or am q can safely be neglected, as we work with a/i q < 0.01 
and om q one order of magnitude smaller than a^ q . Concerning the uncertain- 
ties on ca, cy and or, their statistical errors represent very tiny effects, whereas 
the intrinsic O(a) ambiguity of these coefficients by definition affects any scaling 
observables only at 0(a 2 ). Nevertheless, in the case of cy, for which the avail- 
able non-perturbative estimates at low values of [3 are non-small (of order 0.1) 
and significantly different from each other, one might want to checkQ the effect 
of employing for our observables the one-loop value [20] or the non-perturbative 



estimate by Ref. [13] rather than the value determined in Ref. [19|. 

By definition, among our observables, only myL, r/v and rjy depend on cy. The 
dependence on cy is due to the contribution from [fey(xo)] R to the correlator 
[/cy 1 ,(xo)]R. Since this term comes with a factor of cos a, see eq. (|3.11| ), one can 
expect that for our choice of renormalized quark mass parameters, eq. fl3.2|) , the 



cy-dependent contribution to [ky, (xo)]r is quite small. Indeed, at j3 = 6 we find 
myL = 2.697(18) , r]y = 0.1547(21) , ?jy = 0.1845(21) (4.2) 

with the one-loop value of cy and 

myL = 2.696(18) , r]y = 0.1541(21) , ?jy = 0.1835(21) (4.3) 

with the non-perturbative value of cy given in Ref. [O]. Comparing with the 
values quoted in Table [2] (set A), we see that on myL the effect of using the value 



of cy by Ref. [jl3[ is negligibly small, while the analogous effect on r/y and rjy is 
less than one standard deviation. Moreover, employing the one-loop value of cy 
at j3 = 6.26 we observe deviations from the results of Table that are smaller by 



3 We thank the referee for suggesting this check. 
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about a factor of two than the corresponding deviations found at /3 = 6. We hence 
conclude that the choice of the non-perturbative definition of cy is not important 
for the scaling behaviour of our observables. 

The available non-perturbative estimates of ca are smaller than the corresponding 
ones of cy by at least a factor of four. Among our scaling observables only fhpsL, 
rjPsL and rjpgL depend on ca, and the dependence vanishes as a — > ir/2. The 
situation is hence similar to the case of cy. Concerning the error associated 
with the use of one- loop values for ct, we just remark that the non-perturbative 
estimate given in Ref. |l3j at f3 = 6 is not much larger than the one-loop value, 
which is about 0.02, and affected by large relative uncertainties. 

In view of these remarks and for the sake of simplicity, in our analysis we have 
neglected all uncertainties on the values of the improvement coefficients. 

• Uncertainties on Lm R 

By extra simulations at the same bare parameters as those of the point A in Ta- 
ble [l], but with values of k such that Lm R ~ 0.010 and Lm R ~ 0.034, we estimate 
the derivatives of our meson observables with respect to Lm R . All derivatives are 
of order 1 and compatible with zero within errors. These estimates are employed 
to (slightly) move the central values of the observables, so that the nominal renor- 
malization condition Lm R = 0.020 is exactly matched, and to add quadratically 
to the statistical errors the uncertainties arising from the quoted error on Lm R . 
The effect of this correction is however very tiny, as well as the modification of 
the errors on the meson observables. 

• Uncertainties on L/i R 

By extra simulations at the same bare parameters as those of the point A in Ta- 
ble [l], but with values of fi q such that L/i R ~ 0.140 and L/j, r ~ 0.168, we estimate 
the derivatives of our meson observables with respect to L/i R . All derivatives 
take values between 0.5 and 3, with relative errors less than 10%. The estimated 
uncertainties are employed to add quadratically to the statistical errors the un- 
certainties arising from the quoted error on L/i R . The corresponding increase of 
the statistical errors on the meson observables is significant only in the case of r/ps 
and r^ps- 

• Uncertainties on L/tq 

By an extra simulation at the same bare parameters as those of the point A in 
Table |l|, but for (3 = 6.06 and a value of k such to maintain Lm R <C Lfx n , we finally 
estimate the derivatives of our meson observables with respect to L/tq. We obtain 
estimates of order 1 for the would-be meson masses, of about 0.3 for ?7pg and rjips 
and smaller than 0.1 for rjy and rjy. The relative errors on these estimates are 
of about 10%, except for the derivatives of rfy and rjy which have much larger 
relative errors. Also in this case, the estimated derivatives (or conservative upper 
bounds on them) are employed to add quadratically to the statistical errors the 
uncertainties arising from the quoted error on L/tq. The corresponding increase 
of the statistical errors on the meson observables is typically not larger than half 
the standard deviations in Table g. 
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Set 


'"'Fb- 1 -' 


Uly L 




'IV 


A 


1 sfii n si 

J- .UU1 I 101 




D 5434(53"l 


1 534(29", 


B 


1 8D5('22 N I 


2 652(291 


5570(681 


1565(38') 


C 


1 831(221 


2 6461' 30) 


5514(671 


1615(401 


D 


1.828(29) 


2.652(37) 


0.5499(82) 


0.1596(50) 


set 


m PS L 


fhyL 




rjy 


A 


1.704(12) 


2.395(18) 


0.5686(47) 


0.1823(30) 


B 


1.667(13) 


2.337(22) 


0.5793(53) 


0.1889(39) 


C 


1.680(13) 


2.345(22) 


0.5751(52) 


0.1933(40) 


D 


1.665(16) 


2.324(29) 


0.5756(62) 


0.1944(51) 



Table 3: Final results for our meson observables: the standard deviations given in 
parentheses account for all the uncertainties that we discuss in the text. 

We present the outcome of this analysis in Table |3| and plot the same results versus 
(a/L) 2 in the Fi gures g-|5j. In view of the almost complete implementation of non- 
perturbative 0(a) improvement, for each observable we perform a least squares fit of 
the data to a polynomial of first order in (a/L) 2 . The fit line and the extrapolated 
continuum values with their uncertainty are also shown in the mentioned plots. The 
values of the x 2 per degree of freedom are of order 1 for all the fits we did, and in 
most cases smaller than 1 . The observed scaling violations of our meson observables are 
hence apparently compatible with an 0(a) improved approach to the continuum limit. 

In Table § we quote the extrapolated continuum limit values for all our meson ob- 
servables, together with the relative deviations of these values from the values measured 
at a ~ 0.1 fm (point A of Table |l]). These relative deviations, which can be considered 
as a measure of the size of the lattice artifacts, are indeed fairly small and of the same 
order of magnitude as the analogous relative deviations observed under similar condi- 
tions at a = H|. We remark that the largest relative scaling violation is observed for 
rjy, an observable which was not considered in Ref. || and might be affected by residual 
O(a) cutoff effects due to the poor knowledge of ct- Indeed, both rjy and fhyL, which 
are the only observables that depend on ct, show relative scaling violations larger than 
their cx-independent counterparts, rjy and myL. 



m PS L 


myL 


VPsL 


rjy 


1.801(28) 


2.624(37) 


0.5572(83) [15] 


0.1641(49) [6] 


3.3% 


2.4% 


2.5% 


6.5% 


fhpsL 


fhyL 


rjpsL 


rjy 


1.651(17) 


2.295(28) 


0.5805(65)[15] 


0.2001(50) [7] 


3.2% 


4.4% 


2.1% 


8.9% 



Table 4: Continuum limit values of our meson observables and their relative deviations 
from the values at /3 = 6. The additional errors due to small systematic uncertainties 
in the non-perturbative estimates of Zy and are shown in square brackets. 
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Figure 2: Scaling behaviour and continuum extrapolation for m-p^L and mpgL. 
4.4 Residual O(a) cutoff effects 

The fact that the dependence on a/L of our scaling observables, for the considered 
values of (5 and within relative statistical errors of 1-2%, can consistently be described 
as purely quadratic does not imply the complete absence of cutoff effects linear in a/L. 
We can only conclude that they are small enough to be not clearly visible in our data, 
which in some cases might also be due to accidental cancellations between different O(a) 
cutoff effects. These remarks do not apply to the case of rpcvc> where the cutoff effects 
linear in a/L are expected to be fully negligible, and are actually not seen in the data 
with statistical errors of a few permille. 

On the remaining scaling observables, the only quantitatively significant O(a) ef- 
fects may arise from possibly inappropriate values given either to c± and c\ or to by and 
6a- In the latter case the effects are proportional to a/j, q , which is never larger than 
0.01, whereas the effects arising from the use of perturbative values for a and c\ are 
suppressed in the limit of large T (and are thereby irrelevant for physical applications) . 
The lacking knowledge of the non-perturbative values of these improvement coefficients 
makes any estimate of these uncancelled O(a) effects rather subjective. However, in 
order to disentangle the various residual O(a) effects and get a rough idea of their mag- 
nitude, we have looked at the influence of independent variations of by, b\, ct and c*t 
on our observables. 

When varying by from its one- loop value, eq. ( f4.1|) , to a value of order 1, we observe 
at = 6 a change of a few standard deviations in mps, r?ps and t/ps, and no changes 
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Figure 3: Scaling behaviour and continuum extrapolation for myL and fhyL. 



elsewhere. If we instead vary 6 a from its one-loop value, eq. ( |4.1| ), to a value of order 
1, we note at j3 = 6 a change of about one standard deviation in Lm K , rjy and rjy, and 
negligible changes in all other quantities. 

Varying the values of ct and c\ requires to perform extra simulations. We have 
hence repeated the simulations at the bare parameters of point A and point C of Ta- 
ble |l] by changing either the value of Ct or the value of c\ . Following Ref. || , we have 
chosen the new values of Ct and c\ so that Ct — 1 and c\ — 1 are about 2 and 10 times, 
respectively, larger than the values employed in our scaling test. A variation of ct by 
this amount induces no statistically significant changes in our meson observables. Under 
the mentioned variation of c\, 



c t = 1 - 0.018^ 



ct 



1 - 0.180^ 



we do see statistically significant changes in a few observables, namely rhpsL, ??ps-^ an d 
rjpsL. At f3 = 6 (point A of Table |l|) the observed changes amount to 4-6% of the 
mean values of these observables, whereas at (5 = 6.26 (point C of Table |l|) the relative 
changes are smaller by a factor of 1.5, i.e. they scale proportionally to a/L as expected. 
If one propagates the changes observed at (3 = 6 and f3 = 6.26 to the other values of 
considered in this work (by assuming that they are proportional to a/L), and then tries 
to fit the resulting values of fhpsL, rjpsL and rjpsL for c\ = 1 — 0.180gQ to a constant 
plus a term oc (a/L) 2 , still reasonably good fits are obtained. 

These quantitative checks about the influence of reasonable changes of by, c t 
and ct on the meson scaling observables, together with the possibility of accidental 
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Figure 4: Scaling behaviour and continuum extrapolation for n-p^L and rjpsL- 

cancellations, show that our results in Subsection 4.3 are compatible with the presence 
of residual O(a) cutoff effects, which may individually have a relative magnitude of a 
few percents. On the other end, fitting the results of Table |3| to a polynomial of second 
order in a/L yields continuum limit values that are consistent with those in Table |], 
but with much larger uncertainties. As it is clear from Figures §-|5|, the coefficient of 
the term oc a/L always takes values that are consistent with zero within errors. 

5 Conclusions 

We have presented a scaling test for some representative hadronic observables in the 
pseudoscalar and vector meson channels computed in quenched lattice twisted mass 
QCD with Schrodinger functional boundary conditions. We have also studied an ob- 
servable that allows to quantify the lattice cutoff effects in the PCVC relation, which is 
non-trivial due to /z q ^ 0. To get accurate results at moderate computational effort, an 
intermediate-volume system of physical size 0.75 3 x 1.5 fm 4 has been considered. 

We find that in the parameter region specified by (3 > 6, L/^r = 0.153 ^> Ltur = 
0.020 and T = 2L ~ 1.5 fm the O(a) improvement programme of tmQCD can be 
successfully implemented, although non-perturbative estimates of by and 6a are still 
missing. The size of the observed scaling violations, which range from 0.5% to 9% 
depending on the observables, is acceptably small and comparable to the size of the 
cutoff effects observed in standard lattice QCD with Wilson quarks. Studies of lattice 
tmQCD in larger spatial volume, L = 1.5-2.2 fm, confirm that the pseudoscalar and 
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Figure 5: Scaling behaviour and continuum extrapolation for rjy and rjy. 



vector meson channels can be studied well in the chiral regime by simulations at (3 > 6 



Since this work represents the first non-perturbative study of lattice tmQCD, it is 
worthwhile to emphasize that we employ completely standard lattice techniques and find 
the computational effort to be essentially the same as for Wilson quarks with standard 
mass parameterization. These findings, which have little to do with our particular 
choice of boundary conditions, are basically due to the simple flavour structure of the 
considered correlation functions (see Subsection 2.2) and the good performance of the 
BiCGStab solver in the explored parameter region (see Table |l]). Following Ref. Q, the 
renormalization of the twisted mass parameter turns out to be easy in practice, because 
it can be traced back to the renormalization of the non-singlet pseudoscalar density in 
the massless theory. 

Of course it would be interesting to perform similar scaling tests of lattice tmQCD 
for different choices of the physical parameters: e.g. at L = 0.75 fm but T > 2L, in 
order to render the O(a) improvement of the SF boundaries unimportant, or at higher 
values of /U R ro, as well as for the lattice theory with two twisted light quarks and one 
heavier standard quark, which is relevant for the study of kaon physics p3|J|. Some 
work in this direction is planned for the near future. 



with small scaling violations [11]. 
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